Gravitating SU(N) Monopoles from Harmonic maps 
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Abstract 

Spherically symmetric solutions of the SU(N) Einstein- Yang-Mills-Higgs system are constructed 

using the harmonic map ansatz. This way the problem reduces to solving a set of ordinary differ- 
On. 

ential equations for the appropriate profile functions. In the SU (2) case, we recover the equations 
studied in great detail previously while in the SU(N) (N > 2) one we find new solutions which 
correspond to monopole-antimonopole configurations. 
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I. INTRODUCTION 

Topological defects [jj are thought to have formed during the phase transitions that took 
place in the early universe. If a symmetry is spontaneously broken (G — > H), the topology 
of the vacuum manifold M. = G/H determines which type of the defect arises during the 
symmmetry breakdown. This way we could get domain walls, cosmic strings, monopoles 
and textures if M. is (respectively) disconnected, has contractible loops, non-contractible 
2-spheres and non-contractible 3-spheres. In 1974, 't Hooft and Polyakov [2] realised that 
the bosonic part of the Georgi-Glashow model, which is essentially an £77(2) Yang-Mills- 
Higgs system, possesses soliton solutions which, due to their topological properties, carry 
a magnetic charge. Since the unbroken U(l) group is associated with the electromagnetic 
field, these solutions are said to be describing "magnetic monopoles". Soon afterwards 
people started looking at various embeddings of £77(2) into higher gauge groups 
The embedding into SU (3) was first studied in and the two possible embeddings which 
correspond to a SU(2), respectively £0(3) subgroup, were investigated. For the latter 
embedding, solutions with magnetic charge ±\/3 p as well as solutions with zero topological 
charge [3] were constructed. A systematic analysis of the solutions in an SU (3) model 
with a non- vanishing potential has been done in jjj. Recently, static monopole solutions 
of the second order SU(N) BPS Yang-Mills-Higgs equations, which are not solutions of 
the first order Bogomonlyi equations, have also been constructed j9j. These spherically 
symmetric solutions may be interpreted as monopole-antimonopole configurations and their 
construction involves the use of harmonic maps into complex projective spaces. 

When an ££7(2) Yang-Mills-Higgs system is minimally coupled to gravity, three different 
types of solutions are possible; namely embedded Reissner-Nordstrom solutions Q|, grav- 
itating monopoles and non-abelian black holes Gravitating monopoles exist only up 
to a maximal value of the gravitational coupling beyond which their Schwarzschild radius 
becomes larger than the radius of the monopole core. At this maximal value of the gravita- 
tional coupling the solutions bifurcate producing a branch of extremal Reissner-Nordstrom 
solutions. Non-abelian black holes can be thought of as black holes situated inside the core 
of a magnetic monopole. Consequently, they exist only in a limited domain of the gravi- 
tational coupling-horizon plane. Note that gravitating monopoles in SU(3) corresponding 
to a £77(2) subgroup have been studied in |l2j, while the gravitating monopoles and non- 
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abelian black holes in 577(5), corresponding to a SU(2) subgroup, were studied in jlHI ]. 
Here we will construct the spherical symmetric gravitating solutions of the second order 
SU(N) Einstein- Yang- Mills- Higgs equations which are neither solutions of the first order 
Bogomolnyi equations nor simple embedding of the 577(2) ones. 

When one considers fields in three dimensional space, it is sometimes convenient to intro- 
duce spherical polar coordinates to describe all points in space with the origin of coordinate 
system located at a specific point - like the centre of the soliton (e.g. monopole). Then one 
introduces a radial variable r and two angular variables describing points on the sphere of 
radius r - the complex variables z and z (discussed later in more detail). Therefore, for a 
given r, we have maps from S 2 of radius r, which we can interpret as compactified R 2 . The 
harmonic map ansatz exploits this property; it uses maps of S 2 — > G to construct maps 
of R 3 —>■ G. This is done by considering maps of R 2 — > G and assuming that the parameters 
of these maps are functions of r only. This cannot be done in an arbitrary way; the fact 
that the resultant maps are those of R? — > G imposes some constraints of continuity etc 
(discussed in the next section). 

Our paper is organized as follows. In Section II we present the model and introduce the 
harmonic map ansatz. In Section III we derive and then discuss the resulting equations for 
the SU (2) case while in Section IV the equations for the SU(3) model and also present our 
numerical results. Finally, in Section V we presents our conclusions. 

II. THE MODEL 



The SU(N) Einstein- Yang-Mills-Higgs action is given by: 



S 



R 1 „,,„s 1 T T x 1, / / T o 2\\ 2 



-tr (F^ F>") - -tr (D^ D»$) + -A (tr ($ 2 - V 2 )) ' 



'— g d 4 x (1) 



16ttG 2 v ' 4 

where g denotes the determinant of the metric while the field strength tensor is defined by: 

= d^A u - d v A„ + [An, A u ] (2) 

and the covariant derivative of the Higgs field reads: 

D^ = d^+{A^§}. (3) 

The matrix rj represents a constant matrix of the form: t] = ivl]y, where v G R and ljv 
denotes the unit matrix in iV dimensions. The constants in the action represent Newton's 



constant G, the Higgs self-coupling constant A and the vacuum expectation value of the 
Higgs field v. Note that for N > 2 the potential in (0) is not the most general one that 
could have been used. 

Variation of the action ([TJ with respect to the metric g^ v leads to the Einstein equations 

R^-^vR = 8-kGT^ (4) 
with the stress-energy tensor = g^ v C — 2t^j given by 

+i^A(tr(4. 2 -„ 2 )) 2 . (5) 

Variation with respect to the gauge fields and the Higgs field $ leads to the matter 
equations 

' D^^g-F^-l&D"®} = 0, 



^=g - , 4 

' D M (v/=££>"$) + Atr ($ 2 - rf) $ = . (6) 



In what follows we consider the static Einstein- Yang-Mills-Higgs equations in order to 
construct their spherically symmetric and purely magnetic (ie A = 0) solutions based on 
the harmonic map ansatz first introduced in 



A. The Harmonic Map Ansatz 

The starting point of our investigation is the introduction of the coordinates r, z, z on 
R 3 . In terms of the usual spherical coordinates r, 9, <p the Riemann sphere variable z is given 
by z = e 1 ^ tan(6 l /2) and z is the complex conjugate of z. In this system of coordinates the 
Schwarzschild-like metric reads: 

ds 2 = -A*(r)B(r)dt 2 + -±-dr 2 + 4 f dzdz, B{r) = 1 - (7) 

B(r) (1 + |z| ) r 

where A and B are the metric functions which are real and depend only on the radial 

coordinate r, and m(r) is the mass function. The (dimensionfull) mass of the solution is 

given by = m(oo). For this metric the square-root of the determinant takes the simple 

form: 

_ . . . 2r 2 



Using J7J) the matter equations (jHJ) read: 

i [Z> r $, $] - (1 ^f )2 p 2 F r2 - + D- z F rz ) = 0, (9) 
^ <S>] + ^D r (ABF rz )- 1 ^D z ((l + \z\ 2 ) 2 F zz ) = 0, (10) 

i d 4- id 2 i 2 

£> r (ABr 2 D r $) + 1 \ ] 1 {D Z D- Z <& + D- Z D Z <&) + Atr ($ 2 - ?? 2 ) $ = (11) 



Ar 2 T v 2r 
In addition, the Einstein equations (jlj) take the form: 

im' = SttGT,?, ?^ S = SttG (T ° - 2?) (12) 

where the prime denotes the derivative with respect to r and 

Rfl -L | Z R2 (1 — I — |~|2\4 B 

T ° = - ^ + 1-1) tr [\F rz \ 2 ) + ^^tr (** ) - |tr ((D r *f) 



-tr(|2?/&| a ) + S (tr(* a , (13) 



4r 2 

T ° - T; = - 2E(1 r + 2 N2)2 tr (|F„| 2 ) - |tr ((D r $) 2 ) . (14) 
Next we introduce the following ansatz |9| for the SU(N) gravitating monopoles: 

N-2 , s JV-2 

$ = ^^(P,--), ,1 ^ //; 7> M 7 ; />;. A r = (15) 
i=o ^ ' i=o 

where hj(r), gj(r) are the radial dependended matter profile functions and P(z, z) are N xN 

Hermitian projectors: Pj = Pa = Pf, which are independent of the radius r. Note that all 

N — 1 projectors Pi are orthogonal to each other since P%Pj = for i ^ j and that we are 

working in a realgauge, since A z = —A\. 

As shown in |9J], the projectors defined as 

p *= | Afc/ | 2 ~ . fc = 0,..,/V-l (16) 



A/ = z f - nf ^ ? J) (17) 



where 



give us our required set of orthogonal harmonic maps (for details see |15|). Moreover, the 
harmonic maps with spherical symmetry can obtained by applying the orthogonalization 
procedure to the initial holomorphic vector / given by 



/ (fa fj fx :)'• where fi = J j ^ ( 18 ) 



and ^J 1 ) denote the binomial coefficients. Then equation © is automatically satisfied. 

In dealing with the equations which arise from the harmonic map ansatz (j!5j) it is con- 
venient to replace the profile functions hj(r), gj(r) by the functions bj(r), Cj-(r) which are 
defined as the following linear combinations of g<j and bk'- 

N-2 

hj = ^2 fofc ' c i = 1 - 9j ~ 9j+i, j = 0, . . . , N - 2, (19) 

k=j 

where gN-i = 0. 

Next we will describe in detail the gravitating monopoles obtained from our harmonic 
map ansatz for the simplest cases of SU{2) and SU(3). The situation for general SU(N) 
will then become clear. 

III. GRAVITATING MONOPOLES IN SU(2) 

For N = 2 there are two profile functions, &o,co and our ansatz (JT5j) reduces equations 
(fTU|) and (JTJJ) to the following set of second order nonlinear ordinary differential equations: 

j(ABc' )' = ~b 2 oCo + 1 co (4-1), (20) 
-^(r*ABb' o y = ^b ct + h (bt-^). (21) 

Finally the Einstein equations (|T2*j) take the form: 




where m(r) is given by (J7J). These equations have previously been studied in great detail 
in Q| (after the rescale of the Higgs profile function bo — > 2bo)- We will not repeat the 
numerical calculations here, but refer the reader to the mentioned papers. 

IV. GRAVITATING MONOPOLES IN SU(3) 

For N = 3 there are four profile functions, 6 0) ^i) c 0) c i an d our ansatz (fT5j) reduces 
equations (fTUJl and (jllj) to the following set of second order nonlinear ordinary differential 
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equations: 



-L (r 2 ABi$)' = 1 (26,c| - 6 fc c 2 ) + | b (b) + fcfc, + 6; 



(24) 
(25) 



Here the indices are chosen from the set {0, 1}, k ^ j, and we assume the symmetry under 
the interchange of indices 1 when applied to both the bj and Cj functions. Note that in 
the flat limit, ie for A = B = 1, and for A = the equations of are recovered. 
Finally, the Einstein equations (j!2)) take the form: 



-m 



8ttG 



L f (b'o 2 + W + b? ) + ^ (& 2 C 2 + 6 ?c?) + 2 J (c 2 + c ; 2 ) 



2 A / 

+ ^ ( X - c o - c i + c o + c ? - c o c i) + 7£ ( 6 o + boh + b : 



l--v 2 
1 2 



2 A' 
r A 



8ttG 



\ (b- + M + 6?) + ± (c' 2 + c?) 



(26) 
(27) 



These equations correspond to the (ti, r 2 , r 3 ) = (A 7 , — A 5 , A 2 ) embedding of SU{2) into SU (3) 
(see section below). The Tj's here are the SU{2) (or 5*0(3)) generators, while the Aj's denote 
the Gell-Mann matrices. 

The equations corresponding to the (ri,r 2 ,r 3 ) = (Ai/2, A 2 /2, A3/2) embedding of SU{2) 
into SU(3) have been studied in and are not included in our approach as can be seen 
by setting either c or c\ equal to zero and comparing our equations with those of [3| . This 
is due to the fact that in our construction the corresponding Higgs and gauge fields are 
non-embeddings of the SU (2) ones. 



A. Comparison with the spherically symmetric ansatz 



$ = F 1 (r)Y + F 2 (r) ( Y 2 



The spherically symmetric ansatz used in |7J, |8j, |10( is given as follows: 

2" 



where Y = f • A, A = (A7, — A5, A 2 ) (or in components Y a b = —i^abc^c) and 



, G(r)_ H(r) 



2r v ' 2r 

where [ , ]+ denotes the anti-commutator 



(r x A)j, f • A 



A = 0, 



(28) 



(29) 
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After some algebra, it can be seen that the above Higgs field can be rewritten in terms of 
hermitian and orthogonal projectors constructed out of the matrix Y (note that Y 3 = Y): 

$ = (b + h) (p - + h (Pi - , (30) 

with 

Y + Y 2 

Po = —^, P 1= Y 2 -1. (31) 

With these projectors, the gauge fields can be constructed in an analogous way. Finally, the 
full correspondence between the harmonic map ansatz and the spherically symmetric ansatz 
can be established by observing that 

F 1 = ^(b + b 1 ), F 2 = ^{b -b 1 ), G = G-2 = c - Cl , H = c + Cl . (32) 
B. Numerical results 



As in 



without any loss of generality, we set the vacuum expectation value v = 1 
(different v's can always be accommodated by rescaling the radial coordinate and the Higgs 
field functions) and we also define a 2 = AttG. Then the dimensionless mass of the solutions 
is defined as: M = ^f. Here we will consider only the A = case and leave A ^ to a 
further study. 

In this section we construct the monopole solutions of equations (j24)l - (j27j) numerically 
using a collocation method for the boundary- value ordinary differential equations In 
this procedure the set of non-linear coupled differential equations is solved using a damped 
Newton method of quasi-linearization. 

The boundary conditions for the metric profiles read: 

m(r = 0) = , A(r = oo) = 1 (33) 

and we also assume that A(r = 0) is finite. 

Since the equations (|2*H) - (|2*7j) are very similar to those of the SU(2) case, we expect that 
their solutions will bifurcate into extremal Reissner-Nordstrom solutions at some critical 
value of the gravitational coupling constant a. Prior to describing our numerical results, 
it is worth mentioning that the charge of the limiting Reissner-Nordstrom solution can be 
given in terms of the asymptotic values of the gauge fields c = c (oo), C\ = Cj(oo). Due 
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to (|26|) and taking into account the asymptotic behaviour of the various fields, we find that 
the mass function miyv(r) of the Reissner- Nordstrom solution is given by: 

m RN (r) = nioofiN — , Q 2 = 4 (l - Cq - c\ + 5q + c\ - <^c?) . (34) 

Note that for the extremal Reissner-Nordstrom solution m 00) ^jv — a Q- 

Following the discussions in we note that there are three different types of solutions 
which seem to be of particular interest. They can be distinguished from each other by 
whether they satisfy, or not, the Bogomolnyi equation and whether the symmetry-breaking 
(SB) is maximal (unbroken group U(l) x U(l)) or minimal (unbroken group U(2)). 

First, we will discuss in detail the case of non-Bogomolnyi maximal SB solutions which 
satisfy the following boundary conditions for the fields: 

6 (r = 0)=0, 6 1 (r = 0) = 0, c (r = 0) = 1, Cl (r = 0) = 1 (35) 

at the origin and 

b Q (r = oo) = —2, b\(r = oo) = 4, co(r = oo) = 0, C\(r = oo) = (36) 

at infinity. Any of such solutions has a magnetic charge (0, 2) 9] and according to the 
discussion given there can be interpreted as a superposition of two monopoles and two pairs 
of monopoles-antimonopoles. In the flat limit, this solution has a mass M = 4.5. When 
gravity is minimally coupled to the system, the solution gets smoothly deformed by it. The 
metric function B(r) develops a minimum which gets deeper as a increases. At a critical 
value of a the solution develops a double zero at some finite value of r = rn, which can be 
interpreted as the horizon of the extremal Reissner-Nordstrom solution. Thus, the limiting 
solution can be described by this field for r > r^, while it is non-trivial and non-singular for 
< r < r h . This is illustrated in Fig. 1 (respectively 2), where we present the profiles of 
the metric and gauge functions (respectively, metric and Higgs functions) for the flat limit 
a = 0.0 and a ~ a cr = 0.615. 

As in the £77(2) case, the numerical results indicate that the gravitating solution exists 
up to a maximal value of a, a max ~ 0.625. By decreasing a from a max , another branch 
of solutions can be constructed and B min = is reached at the critical value a cr ~ 0.615. 
This is illustrated in Fig. 3, where plots of and of the minimal value B m i n in terms 
of a are presented. Note that at a = a cr , the quantity m ^ L equals two. Since for the 
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Reissner-Nordstrom solution m °° Q flAr = Q, this implies that the solution indeed bifurcates 
into a charge-two Reissner-Nordstrom solution (as can be calculated from (1341) ). 

Q 

Next we studied the gravitating analogue of the Burzlaff solution p|. In this case bo = —b\ 
while Co = Ci and the solution corresponds to a non-Bogomolnyi non-maximal SB case with 
charge (0, [2]). [Here the notation is that magnetic weights are defined by square brackets 
We choose the following boundary conditions for the fields: 

b (r = 0) = 0, &i(r = 0) = 0, c (r = 0) = 1, c x (r = 0) = 1 (37) 

at the origin and 

bo{r = oo) = \/3, b\{r = oo) = — v3, co(r = oo) = 0, ci(r = oo) = (38) 

at infinity. The choice of the boundary conditions for the Higgs field at infinity is not 
fixed for a vanishing potential; however, following we have chosen our conditions as if 
the potential was present. This solution, due to (JM)), should bifurcate producing a charge- 
two Reissner-Nordstrom solution as confirmed in Fig. 3. Again, we find a back-bending 
since the gravitating Burzlaff solution exists up to a max ~ 1.106 and reaches B m i n = at 
a = a cr ~ 1.063. In the flat limit, the mass of the solution is M = 2.5. Since the solutions 
are less heavy than the ones in the non-Bogomolnyi maximal SB case, it is clear that the 
critical value of the gravitational coupling should be bigger in the Burzlaff case (which is 
confirmed by our numerical results). 

Finally, we constructed the gravitating analogues of the Bais solution 6] by choosing the 



boundary conditions (again) as if the potential were present [8j. The boundary conditions 



were 



6 (r = 0) = 0, 6i(r = 0)=0, c (r = 0) = 1, a{r = 0) = 1 (39) 
at the origin and 

b (r = oo) = — \/3, b\(r = oo) = 0, c (r = oo) = 0, ci(r = oo) = — -=. (40) 

v2 

at infinity. These solutions correspond to the Bogomolnyi non-maximal SB solutions with 
charge (2, [1]). Due to (|3*4*j) the branch of solutions should bifurcate producing a branch of 
Reissner-Nordstrom solutions with charge v3 as shown in Fig. 3. In fact the gravitating 
Bais solutions exist up to a maximal value of a: a max ~ 1.446. We were not able to find 
a back-bending and it is likely that there is none since B min w 4.7 ■ 10~ 3 at a = a max . If 
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there were a back-bending, then the branch of solutions would be very small. Since, in the 
flat limit, the mass of the gravitating Bais solution (again) is smaller than in the other two 
cases studied, the maximal value of a for the former is larger than for the latter. 



V. CONCLUSIONS 

We have studied gravitating SU (N) monopoles by relying on the harmonic map ansatz. 
In the SU(2) case, we have recovered the equations studied previously in great detail in 
while in the SU(3) one, we have found the gravitating analogues of the solutions obtained by 
the embedding of SO (3) into SU(3). Since in the case of the vanishing potential (considered 
here), the boundary conditions of the Higgs profile functions at infinity are not fixed, they 
can be chosen at will. In fact they have been chosen so that the gravitating analogues of 
the solutions discussed in [9] can be constructed. These solutions correspond to monopole- 
antimonopole configurations and not to single monopoles as constructed in |12| . 

In all three cases studied here, it has been found that the solutions bifurcate producing a 
branch of extremal Reissner- Nordstrom ones with charge Q, which is fixed by the asymptotic 
values of the gauge fields. Interestingly, our numerical results indicate that a second branch 
of solutions, which extends backwards from the maximal possible value of the gravitational 
coupling, exists only in the case of non-maximal symmetry breaking. In the case of maximal 
symmetry breaking (even in the limit of vanishing Higgs coupling) no second branch was 
found. This is in contrast with the case of single SU (2) 11( and SU(3) 12j monopoles where 
the second branch of solutions exists for vanishing or small Higgs coupling. 

In this paper we have not constructed solutions for non- vanishing potential or non-abelian 
black hole solutions. A systematic study of these solutions is left to a future work [18]. In 
particular, we would like to study the corresponding configurations for intermediate Higgs 
coupling constants. It is known that for the single SU(2) and SU(3) monopoles the so- 
called "Lue- Weinberg" |l9j phenomenon was observed; ie for intermediate values of the 
Higgs coupling the solutio ns develop a second "inner" horizon and in the limit of critical 
gravitational coupling describe "hairy" black holes. It will be interesting to see whether this 
phenomenon persists for monopole-antimonopole configurations. 
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FIG. 1: The profiles of the metric A(r) and B(r) as well as of the gauge fields G(r) = co(r)— c±(r)+2 
and H{r) = co(r) + c±(r) are presented for the non-Bogomolnyi maximal SB solutions in the flat 
limit a = 0.0 (dashed) and close to the critical limit (solid). Note that for a = 0.0, A(r) 

and B(r) are constant since A(r) = B(r) = 1. 
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FIG. 2: Same as Fig. 1 for the profile functions of the Higgs field F\(r) = ^(^oM + ^i( r )) an d 
F2(r) = ^(bo(r) — b\(r)). For comparison, we show again the metric functions A(r) and B(r). 
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FIG. 3: The quantity aM = and the minimal value of the metric function B(r), B m i n , in 
terms of a is plotted for a) the Bogomolnyi non- maximal SB case (gravitating Bais solution) , b) the 
non-Bogomolnyi non-maximal SB case (gravitating Burzlaff solution) and c) the non-Bogomolnyi 
maximal SB case. 
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